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Abstract 

The directed percolation (DP) hypothesis for stochastic, range-4 cellular au- 
tomata with acceptance rule y < Ylj=-A s i-j — 6, in cases of y < 6 was 
investigated in one and two dimensions. Simulations, mean-field renormal- 
ization group and coherent anomaly calculations show that in one dimension 
the phase transitions for y < 4 are continuous and belong to the DP class for 
y = 4, 5 they are discontinuous. The same rules in two dimensions for y = 1 
show 2 + 1 dimensional DP universality; but in cases of y > 1 the transitions 
become first order. 
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I. INTRODUCTION 



The phase transition theory of nonequilibrium systems is an area of considerable interest. 
There is not such a clear view about critical universality classes as there is for equilibrium 
systems. The cellular automata (CA) approach for modelling nonequilibrium systems can 
be regarded as a general tool like differential equations. It introduces a space and time 
discretization — thus making the solution treatable — but preserving the features of interest. 
As was shown in |] CA models can describe very rich behavior. Even the simplest one 
dimensional CA can have various phase transitions, can generate chaotic patterns, or show 
complex self-organizing or reproductive behavior. 

The directed percolation (DP) hypothesis @ states that all continuous transitions of a 
scalar order parameter to a unique absorbing state in dimension d belong to the DP uni- 
versality class, represented by the d + 1 dimensional directed percolation || or reggeon field 
theory. There are some additional conditions - usually not stated explicitly - to be fulfilled 
as well : short range interaction both in space and time, translational invariance, absence 
of multicritical points, and the non-vanishing probability for any active state to die locally. 
Although there is no proof for this conjecture there exists a great number of models, which 
have been found to belong to this class. Attempts to find counterexamples |§-|B| by other 
simulation results [|7| have been unsuccessful. The motivation of our investigation is to find a 
possible exception to the DP conjecture in the cellular automata phase transition of different 
dimensions. The stochastic cellular automata - investigated by Bidaux, Boccara and Chate 
(hereafter referred to as BBC CA) || and Jensen using simulations - in one dimension 
showed controversial results concerning the universality class. Steady state simulation gave 
non-DP whereas the time-dependent simulation resulted in the DP class. Even though the 
latter method is considered to be more precise, arguments based on universality suggested 
that the same rule in two dimensions should produce a transition of 2+1 dimensional DP 
class. Simulations proved this to be a first order transition ||. In order to clarify the 
situation we used various selection of different analytical and simulation methods. 

In addition we were interested in the change of universality by varying the lower accep- 
tence value y of the rule : 



By setting up the mean-field steady state equation for the probability of ones (c) : 



the leading term in the polinomial function /, which determines the critical behavior as 
p —>■ p c (and c — ► 0) is c for y = 1, c 2 for y = 2 and c 3 for y = 3. In analogy with the 
equilibrium Ising model — where the three-spin Ising model has a different universality from 
the two-spin nearest neighbor Ising model — we expected deviations from the DP behavior. 
To check this assumption we extended our investigation to the y = 1, 2, 3 cases in one and 
two dimensions. 




where X e {0, 1} is a two valued random variable such that 



Prob(X = 1) = p 



dc 
dt 



f(p,c,t) = 



(1) 



2 



II. STEADY STATE SIMULATIONS 



The steady state Monte Carlo simulations (STDMC) were carried out on L = 40000 and 
L = 256 x 256 sized lattices with periodic boundary conditions. In each case measurements 
at 30-50 different p values were taken. Following 10 5 - 5 • 10 5 relaxation lattice updates 
at each p the concentration averages and r.m.s. values were determined. The exponent (3 
and p c were fitted together - using the weighted least squares error method - to the critical 
scaling formula 

c(p) ~ \p~Pcf ■ (2) 

Since the concentration in the immediate neighborhood of the critical point has a rela- 
tively large statistical error - owing to the increased relaxation time - data points within the 
1% neighborhood of p c were discarded. In order to compensate this truncation of data we 
performed the fitting in several different neighborhoods of the p c . Since we don't know how 
the (|2|) scaling breaks down away from p c , we have checked this tendency of the fitting by 
varying p cut of the used data points : (p c x 1.01) < {p} < p cut . This kind of analysis enabled 
us to see how the critical exponent converged as p cut — > p c and we were able to extrapolate 
for the true (3 exponent in spite of the growing statistical errors in the vicinity of p c . 

A. Results in one dimension 

The generally accepted 1+1 dimensional DP critical exponets are summarized in Table 
|. Figure [l] shows that the estimated critical exponent of the order parameter converges to 
the value f3 ~ 0.28, which is about the DP value determined by Dickman and Jensen by 
series expansion ||. The critical point and (3 estimates ( given in Table ||) were determined 
by linear weighted least-square fitting to the functions of p cut . The results of the BBC [|5| 
(y = 3) model are the most uncertain, but still seem to extrapolate to the DP class value 
if Pcut Pc- Estimates based on fitting with less than 10 data points should not be taken 
very seriously and are shown for information only. 

As a tendency the transitions become more and more steep (and p c — > 1) as we increase 
the lower limit of the CA rule (y). For y = 4 and 5 non-empty steady state exists in the 
deterministic, p = 1 limit only. One can easily check that only a few oscillating structures 
can survive. There is a first order transition to the empty state if p < 1, which is quite rare 
among one dimensional CA. For y > 5 every cell dies out. 

To verify our results further simulations and analytical calculations were carried out; 
these are discussed in the subsequent sections. 

B. Results in two dimensions 

In two dimensions for y — 1 the transition seems to be continuous and the estimates 
for the critical exponent converged to the 2 + 1 dimensional DP class value : (3 ~ 0.58 
(Figure ffl/d). For y = 2 and 3 the transitions appear to be discontinuous, see Figures If 
formula (2) is fitted to the concentration data - as if continuous transitions were assumed - 
results in p c -s with values lower than p with non-zero densities. Since it is not easy to prove 
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first order transition by simulation, we applied the generalized mean-field method, which 



has already been found to be very useful for deciding the order of the transition at CA [JTO 
The results of these calculations are also shown in Figure Q and they are discussed in section 

E 

C. Finite-size scaling analysis 

Finite-size scaling is shown to be applicable to continuous transition to an absorbing 



state of nonequilibrium systems [11]. At the critical point the steady state density (c) and 



the fluctuation x — L d (< c 2 > — < c > 2 ) scale with the system size as : 

c(L)<xL-M"±, (3) 

X (L) <x irrt"\ (4) 
where u± is the correlation length exponent in the space direction : 

£(p) a: \p-p c \~ v± - (5) 

Simulations were done in one dimension for lattice sizes : L = 16,32,64, ...,8192 with the 
corresponding timesteps : t = 200, 400, 800, 102400. Averaging was performed on the 
'surviving' samples out of N s = 262144, 131072, 65536, 512 independent runs. The p c 
values are taken from the time-dependent MC calculations. Figure |3| shows that the c and 
X results fall nicely on a straight line on a log- log plot, but for y = 1 the fitting resulted in: 

P/v ± = 0.260(4), 7 /i/j_ = 0.43(1), 

for y = 2 : 

(3/u ± = 0.252(1), j/p ± = 0.548(7), 

and for y = 3 : 

P/v± = 0.258(1), 7/z/ ± = 0.496(4), 

which differ slightly from the DP class values jO : 

P/v± = 0.2522(6), 7/i/i. = 0.496(2). 

We have shown that the deviations come from systems with small sizes, similarly as for 
. By omitting points with sizes less than L = 16,32,64 the results converge to the DP 



values (see Tables III, |V[ 0) and we can see the emergence of a DP like behavior in the 
large L limit. 

III. TIME-DEPENDENT SIMULATIONS 

Time-dependent Monte Carlo simulations (TDMC) have been shown to be an efficient 
method for locating critical points and estimating exponents for models with continuous 



transition to an absorbing state [|i2Hl4|j . The initial state of this simulation is the smallest a 
single live 'seed', that can survive at the center of the lattice. The evolution of the population 
is followed up to several thousands of lattice update steps and the following quantities are 
measured at each time step : 
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• survival probability s(t) 

• concentration c(t) 

• average mean square distance of spreading from the center R 2 (t) 

The evolution runs are averaged over N s independent runs for each different value of p in 
the vicinity of p c ( but for R 2 (t) only over the surviving runs ). At the critical point we 
expect these quantities to behave in accordance with the power law as t — > oo, i.e. 

s(t) oc t~ s , (6) 
c{t) oc f , (7) 

R 2 (t) oc t z . (8) 

We carried out simulations for the one-dimensional systems up to t = 5000 time steps with 
N s ~ 10 6 independent samples. The seed was a single 1 for y — 1, a neighboring pair of 1-s 
for y = 2, and a triplet for y = 3. We did not follow the evolution of R 2 (t). To estimate the 
critical exponents we determined the local slopes : 

_, (t) = Ms(')M'/H] 

m(m) 
ln(m) 

using m = 8. Figures f|, §, |6| show these quantities as functions of t~ l . In the case of power- 
law behavior we should see a straight line as t^ 1 — > 0, when p = p c . The off-critical curves 
should possess curvature. Curves corresponding to p > p c should veer upward, curves with 
p < p c should veer down. The critical exponent should be read off as the interception of the 
critical (straight) curve with the ordinate axis. 

It is emphasized that it is very difficult to get precise estimates on the critical data by 
visual inspection of the curves . While for y = 3 the middle curve - which seems to be the 
nearest to the critical one - one can read off DP exponents, for y = 1, 2 one could conclude 
that we have non-DP values. The logarithmic c(t) and s(t) data were fitted by parabolas : 

ln(c(t)) = const. + f] X ln(t) + r] 2 ln 2 (t) (11) 

ln(s(t)) = const. + 5 X ln(t) + 5 2 ln 2 (t) (12) 

at different p values. Estimations for p c are based on the conditions 1] 2 = and 5 2 = 0; the 
results are shown in Table |I|. It should be noted, however that even though the difference 
is small (Ap < 10 -4 ) the conditions rj 2 = and 5 2 = do not give exactly the same p c . Also 
the estimates for the exponents are slightly off the DP class values. The errors must have 
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statistical nature with finite size errors being excluded since the fronts of the populations 

— over the entire run — cannot exceed the boundary of the lattice. 

To get more precise exponents we used the general form - including corrections to scaling 

- of the local slopes |H| 

5(t) = 5 + at' 1 + b6T 5 ' + ... (13) 

where 5' is a correction to the scaling exponent. As indicated in Table [II] the exponents 
agree well with the DP class data. 



IV. MEAN-FIELD RENORMALIZATION GROUP ANALYSIS 

Renormalization group methods for nonequilibrium systems have been successful when 
combined with mean-field approximation. At the critical point the scaling behavior of quan- 
tities calculated on small clusters of different sizes is exploited. One main branch of such 
calculations is the Coherent Anomaly Method (CAM) its application to these CA models will 
be discussed in section [V|. Another important approach is the so-called mean-field renormal- 



ization group (MFRG) method, originally introduced by Indekeu et al. |L5| for equilibrium 
systems. Later it was shown that the method may also be applied to self-avoiding walk and 
percolation problems [[l|J. The application to nonequilibrium systems was demonstrated by 
fT7|l . Although this method is not considered to be the most accurate, it is relatively simple 
and effective. Developments to enhance its performance are under way |1| . Here, we utilize 
the approach for one- dimensional CA. 

First we set up mean-field equations for different sized block probabilities (Pi) in the 
steady state. Here the index i denotes an n block. However we do not use the Bayesian 
extension trick as in the case of the generalized mean-field calculations, but take into account 
the external field with the average concentration (c): 

= ^ = /({Pi(p)},c) (14) 

Because in the vicinity of the second-order phase transition to the absorbing state c — > 0, 
we expand and keep linear terms only in the equations. We express the probability of a site 
being occupied (O n (p, c)) by calculating it in terms of different (n) cluster sizes. For the 
y — 1 model in the n = 1 case the equations are : 

-8cpP + (1 -8 c) (1 -p) P t = 0, (15) 



Po + Pi= 1. (16) 
In view of this, the probability of a site being in state 1 is : 

1 (p,c) = ^ . (17) 

K ^ J l-8c-p+16cp K J 

Expanding in factors of c we get : 
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O 1 (p,c) = ^ + O(c 2 ). (18) 

For clusters of pairs (n = 2) the equations look like : 

= -2epP 00 - 12c (1-p) pP 00 -Qc'p 2 P o 
+ 2 (l-8c) {l-p) 2 P 01 
+ (l-8c>) (l-p) 2 Pii, 



= 2cpP 00 + 12c {l-p)pP o 
-2 (l-8c) (l-p) 2 P w -2 (l-8c) p 2 P , 
+ 2 (l-8c) (l-pjpPw, 



= 6cp 2 P 00 + 2 (l-8c) p 2 P 0J 
- (l-8c) (l-p) 2 P„ 
-2 (l-8c) (1-p) pPii, 



Poo+2P M +P w = l, (19) 
Now, the probability of 1 after linearization is : 

^ g ) = 7i (7 x + ? + °^ (20) 
(-i +p) 

At the critical point the scaling requires that the O n (p,c) probability must scale like the 
external field : 

O n '(P, c ') = c _ (21 ) 
O n (p,c) c 1 ' 

This condition enables one to combine the expressions for O n and O n ' and to get an estimate 
for the p c . Comparing the n — 1 and n = 2 results we obtained p c (l, 2) = 0.123106. 

The exponent u±, which measures the critical behavior of the correlation length, can be 
obtained from the scaling of (p — p c ) resulting from change of the length scale I = n'/n : 



dO n (p, c) 



dO n (p, c) 



= l 1 ^ (22) 

P=Pc 



For n — 1 and n = 2 clusters it is : u±(l,2) = 5.802. This is very far from the DP 
class u± ~ 1.1 value, but if calculations are done on larger clusters the results converge 
there. We were able to go up to n = 5 clusters with the help of symbolic MATHEMATICA 
programming. Table shows that there is a slow convergence of the p c and u± values. 
With increasing the cluster sizes we expect more accurate results. 

Unfortunately we were not able to apply the method for y = 2,3 models because the 
mean-field approximations predict first order transitions for small cluster sizes and thereafter 
no scaling is expected. 
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V. GENERALIZED MEAN-FIELD RESULTS 



The generalized mean-field approximation (GMF) introduced by ]TP| , |2T| has already been 
applied to describe phase transitions of different stochastic CA models [f2~T| , |10f with success. 
The phase structure and the tricritical point of the one-dimensional BBC model (y = 3) 



with site-exchange mixing has been explored with the help of GMF [10 



The traditional mean-field approximation has been extended to handle clusters of sizes 
n > 1 with the help of the Bayesian extension. Correlations of range > n at a given 
level of approximation (n) are neglected and the convergence of the results is analysed by 
increasing n; the details of the technique are given in PTJ. The extension of the method to 
the two-dimensional CA model is straightforward. In this case the nearest- and next-nearest 
neighbor of the lattice point can contribute to the sum defined in (1). Taking into account 
the probability of a given distribution we used maximal overlap of the clusters. At three- 
point level the form of the cluster breaks the x — y symmetry and prevents the configuration 
from being covered symmetrically. Because of this the result of the equations depends on 
the way of covering the configuration. To avoid this uncertainty we restricted ourselves to 
symmetrical clusters (n = 1, 2, 4). 



A. Results in one dimension 



In the case of the y = 1 model the convergence of the approximations is monotonic. 
Even the traditional (n = 1) mean-field approximation predicts continuous phase transition. 
By applying quadratic fitting to the p c {n) results we got an extrapolation value of n — > 
oo value, which agrees with the simulation result (see Table |V1J| ). For y = 2,3 the low 
level approximations gave a discontinuous phase transition, but the gap sizes decrease with 
increasing n and finally disappear at n = 6 for y — 2, and at n = 3 for y — 3. It is interesting 
that the convergence is non-monotonic for the y = 3 case in contrast to what we expect for 
equilibrium systems. 



B. Results in two dimensions 



As one expects in two dimensions the GMF approximation gives more precise results than 
in one dimension. For the n = 4 level approximation the p c data differ from the simulations 
by a few per cent only. For the y = 1 model all approximations predict a second order phase 
transition, whereas for y — 2, 3 cases the transitions are discontinuous (see Table |V111| ) and 
the gap sizes increase with n. We can expect that in the n — > oo exact limit the transition 
remains first order in agreement with the STDMC simulations. 



VI. COHERENT ANOMALY EXTRAPOLATION 



The GMF approximations give slow convergence at the critical point since the correla- 
tion length there goes to infinity. Therefore the accuracy of method breaks down in the 
immediate vicinity of p c . To obtain estimates for the critical exponents one needs to apply 
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an extrapolation technique. This technique usually is based on finite size scaling theory. 
The mean-field renormalization technique is one such approach. 

The coherent anomaly method (CAM) introduced by Suzuki |22|] is based on the scaling 
relation that the solution for singular quantities at a given in) level of approximation {Q n {p)) 
in the vicinity of the critical point is the product of the classical singular behavior multiplied 
by an anomaly factor a(n), which becomes extremely large as n — > oo (and p™ — > p c ): 

Q n ~a{n){plp n c -iy>* , (23) 

where p is the control parameter and to c \ is the classical critical index. The divergence of 
this anomaly factor scales as 

a(n) ~ (p n c - p e y~** , (24) 

thereby permitting the estimation of the true critical exponent oo, given a set of GMF 
approximation solutions. The method has been tested on a great number of solved and 



unsolved systems The application to cellular automata was shown by [2l|. We followed 
the technique described there in the case of the present CA models. A new parametrization 
is used suggested by ^5 



Sn = (Pc/p n c ) 1/2 -( P :/Pc) 1/2 , (25) 

that has an invariance property: p <-» p~ x . We took into account a correction term to the 
anomaly scaling: 

a(n) = b 5t^ cl + c 5t^ cl+1 ■ (26) 

The GMF calculation gave enough data to do the CAM extrapolations for the one dimen- 
sional y = 1 case (Table |TXD . In a similar way to the MFRG case for y = 2, 3 models the 
first order predictions of the low level GMF results excluded CAM calculations. Table [X| 
shows that (3 results are stable in contrast to neglecting points from the nonlinear fitting 
procedure by Eq . (|26|) and indicating DP class behavior. 



VII. CONCLUSIONS 

A deatiled analysis for a family of totalistic, critical CA in one and two dimensions has 
been carried out in order to check the DP conjecture and to clarify the former contradictory 
simulation results of the BBC (y = 3) model. Different critical exponents were determined 
by simple steady state, finite size scaling and time dependent simulations. The simulation 
results are more extensive than the previous ones and none of them seems to break the 
DP rule. Rather they all give estimates for different critical exponents that are close to 
the DP universality values. We wish to emphasize that careful analysis is necessary for 
each kind of simulation to get unambiguous results. We have shown that the systematic 
calculation of GMF provides a reliable way to distinguish the order of phase transitions. 
Mean-field renormalization group and coherent anomaly method extrapolations to GMF 
provided estimates for critical exponents. These results are again in accordance with the 
DP conjecture. 
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FIGURES 



FIG. 1. Convergence of the (3 exponent estimate in the scaling region. Data are from steady 
state simulations. The number of points corresponds to p cu t such that the distance between sub- 
sequent abscissa points is Ap ~ 0.001. (a) One dimension y = 1. (b) One dimension y = 2. (c) 
One dimensions y = 3. (d) Two dimensions y = 1. 

FIG. 2. Concentration versus p in two dimensions for (a) y = 2.; (b) y = 3 models. Dashed 
lines correspond to GMF results of n = 1,4 point approximation. Data with error bars come from 
steady state simulation 

FIG. 3. Finite size scaling of c(L) and x{L) for the 1 dimensional BBC (y = 3) model. For 
y = 1, 2 similar plots were obtained. 

FIG. 4. Time dependent simulation results in one dimension for the y = 1 model at p values 
(from bottom to top) : 0.22086, 0.2209, 0.22092, 0.22095, 0.22098. (a) The local slope of the 
logarithm of the c(t) versus ln(l/i) is plotted. The intercept of the scaling curve with the ordinate 
gives exponent i]. (b) The local slope of the logarithm of the s(t) versus hi(l/i) is plotted. The 
intercept of the scaling curve with the ordinate gives exponent 5. 

FIG. 5. Time dependent simulation results in one dimension for the y = 2 model at p values 
(from bottom to top): 0.4354, 0.43545, 0.43547, 0.4355, 0.4356. (a) The local slope of the log- 
arithm of the c{t) versus ln(l/t) is plotted. The intercept of the scaling curve with the ordinate 
gives exponent rj. (b) The local slope of the logarithm of the s(t) versus ln(l/i) is plotted. The 
intercept of the scaling curve with the ordinate gives exponent 5. 

FIG. 6. Time dependent simulation results in one dimension for the y = 3 model at p values 
(from bottom to top) : 0.7216, 0.7218, 0.7220. (a) The local slope of the logarithm of the c(t) 
versus ln(l/t) is plotted. The intercept of the scaling curve with the ordinate gives exponent rj. (b) 
The local slope of the logarithm of the s(t) versus ln(l/t) is plotted. The intercept of the scaling 
curve with the ordinate gives exponent 5. 
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TABLES 



TABLE I. DP class critical exponents 



f3 rj 5 u± 

0.2767(4) 0.3137(1) 0.1596(4) 1.0972(6) 



TABLE II. Simulation results in one dimension. 





y = l 


y = 2 


y = 3 


TDMC 








Pc 


0.22095(3) 


0.43543(5) 


0.7217(2) 


V 


0.317(7) 


0.317(7) 


0.307(9) 


5 


0.149(5) 


0.159(7) 


0.161(8) 


STDMC 








Pc 


0.22103(4) 


0.435437(9) 


0.7198(1) 


(3 


0.267(2) 


0.258(1) 


0.298(7) 



TABLE III. Small size effects in FSS as a funtion of cutoff for y = 1 in one dimension 

cutoff 16 32 64 DP 

P/v± 0.260(4) 0.257(1) 0.255(1) 0.251(1) 0.2522(6) 

j/u ± 0.43(1) 0.443(6) 0.462(6) 0.48(1) 0.496(2) 



TABLE IV. Small size effects in FSS as a funtion of cutoff for y = 2 in one dimension 



cutoff 





16 


32 


64 


128 


DP 


(3/v± 


0.252(1) 


0.254(1) 


0.254(2) 


0.253(3) 


0.252(4) 


0.2522(6) 


l/v± 


0.548(7) 


0.537(9) 


0.52(1) 


0.51(1) 


0.50(2) 


0.496(2) 



TABLE V. Small size effects in FSS as the funtion of cutoff for y = 3 in one dimension 



cutoff 





16 


32 


64 


DP 


(3/v ± 


0.258(1) 


0.258(1) 


0.255(1) 


0.253(3) 


0.2522(6) 




0.496(4) 


0.497(4) 


0.493(7) 


0.49(1) 


0.496(2) 
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lAoLiiii VI. 


Mi ixKj results tor y = 


1 in one dimension 


(ni + n i+1 ) 1 




p c (i,i + 1) 




v±(i,i + 1) 


1/3 




0.123106 




5.802 


1/5 




0.135359 




3.287 


1/7 




0.147078 




2.222 


1/9 




0.157571 




1.627 


expected values 




simulation 




DP class 






0.22095(3) 




1.0972(6) 


TABLE VII. 


Convergence of the critical point estimates of the 


one dimensional CA calculated 


by GMF approximation. First order transitions are denoted by boldface numbers. Gap sizes (c(p c ) 


are shown for y - 


= 2,3. 








n 


y = l 


y = 2 




y = 3 




Pc 


Pc 


c(Pc) 


Pc C(p c ) 


1 


0.111 


0.354 


0.216 


0.534 0.372 


2 


0.123 


0.378 


0.228 


0.649 0.385 


3 


0.135 


0.401 


0.236 


0.903 0.0 


4 


0.146 


0.419 


0.215 


0.800 0.0 


5 


0.157 


0.427 


0.096 


0.736 0.0 


6 


0.166 


0.426 


0.0 


0.703 0.0 


extrapolation 


0.196(6) 


0.467(7) 


0.0 


0.72 ±0.1 0.0 


simulation 


0.22095(3) 


0.43543(5) 


0.0 


0.7217(2) 0.0 


TABLE VIII. 


Convergence of the critical point estimates of the two dimensional CA calculated 


by GMF approximation. First order transitions are denoted by boldface numbers. Gap sizes (c(p c ) 


are shown for y = 


= 2,3. 








n 


y = l 


y = 2 




y = 3 




Pc 


Pc 


c(Pc) 


Pc C(p c ) 


1 


0.111 


0.354 


0.216 


0.534 0.372 


2 


0.113 


0.326 


0.240 


0.455 0.400 


4 


0.131 


0.388 


0.244 


0.647 0.410 


simulation 


0.163 


0.404 


0.245 


0.661 0.418 
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TABLE IX. The basis of CAM fitting for y = 1 in one dimension 



n 5 n a(n) 

1 0.701078 0.248755 

2 0.593013 0.331049 

3 0.494564 0.428537 

4 0.409735 0.546581 

5 0.339569 0.673658 

6 0.284204 0.810893 



TABLE X. Stability of (3 exponent for y = 1 in one dimension calculated by CAM 

data set (3 

{1,2,3,4,5,6} 0.273 

{ ,2,3,4,5,6} 0.280 

{1 ,3,4,5,6} 0.273 

{1,2 ,4,5,6} 0.283 

{1,2,3 ,5,6} 0.260 

{1,2,3,4 ,6} 0.273 

{1,2,3,4,5 } 0.266 
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